Rogue decoherence in the formation of a macroscopic atom-molecule superposition 



Olavi DannenbergJ 

Helsinki Institute of Physics, PL 64, FIN-00014 Helsingin yliopisto, Finland and 
Department of Physics, University of Turku, FIN-20014 Turun yliopisto, Finland 

Matt Mackie 

Department of Physics, Temple University, Philadelphia, PA 19122 and 
Department of Physics, University of Connecticut, Storrs, CT 06268 
(Dated: February 2, 2008) 

We theoretically examine two-color photoassociation of a Bose-Einstein condensate, focusing on 
the role of rogue decoherence in the formation of macroscopic atom-molecule superpositions. Rogue 
dissociation occurs when two zero-momentum condensate atoms are photoassociated into a molecule, 
which then dissociates into a pair of atoms of equal-and-opposite momentum, instead of dissociating 
back to the zero-momentum condensate. As a source of decoherence that may damp quantum 
correlations in the condensates, rogue dissociation is an obstacle to the formation of a macroscopic 
atom-molecule superposition. We study rogue decoherence in a setup which, without decoherence, 
yields a macroscopic atom-molecule superposition, and find that the most favorable conditions for 
said superposition are a density p ~ 10 12 cm -3 and temperature T ~ 10~ 10 K. 

PACS numbers: 03.75.Gg, 03.65.Ta, 32.80.Wr 



I. INTRODUCTION 

Since the famous thought experiment by Schrodinger 
, macroscopic superpositions have been a part of quan- 
tum folklore that have annoyed specialists of the field. 
Schrodinger considered a cat that was sealed into a steel 
chamber with a diabolic device that consists of a small 
amount of radioactive material and a Geiger counter that 
was attached to a hammer device. If one radioactive nu- 
cleus would decay, then the Geiger counter would trigger 
the hammer device, and the hammer device would scrap 
a small bottle of Prussic acid which kills the cat. If the 
problem is studied strictly quantum mechanically, from 
outside the steel chamber, one finds that after a while the 
cat would be in a superposition state of alive and dead. 
If the cat would be in the statistical mixture of being 
alive or dead, the situation would not differ much from, 
e.g., tossing a coin with the probability of 0.5 for heads 
or tails. In a classical model, we do not yet know the 
outcome, but the state of the cat (coin) is either alive or 
dead (heads or tails), even before the measurement. The 
quantum cat (coin) behaves differently, with laws that 
allow superpositions between the outcome alive and out- 
come dead (heads and tails). Hence, rather than certain 
alive or certain dead (heads or tails) before the measure- 
ment, there must be quantum correlations between the 
the two outcomes. 

For Schrodinger, the reason for putting the cat in the 
experimental spotlight was obvious. The absurd outcome 
of the experiment illustrated that, if the outcome was a 
result of treating everything-including the cat-strictly 
quantum mechanically, something was wrong, if not in 



the theory itself, then at least in the understanding quan- 
tum mechanics. As previous decoherence studies have 
shown 0, H H H H 0, |H, one should be wary of too 
many oversimplifications and approximations. The en- 
vironment of the cat should also be modelled, not only 
a point-like cat in a vacuum. The interaction between 
the cat and the environment induces decoherence that 
damps quantum correlations out of the cat, resulting in 
a classical statistical mixture of living and dead cat. In 
other words, the cat is too macroscopic to be in a su- 
perposition state. But what, then, is the upper limit of 
macroscopicity for a superposition state? 

Recent studies [1 H GH H 13 13 have pointed 
out that the macroscopicity question can perhaps be 
addressed with Bose-Einstein condensates. References 
Ei EH El consider superpositions of two-component 
atomic condensates, while dynamical creation of a macro- 
scopic superposition between atomic and molecular con- 
densates via photoassociation 0] or Feshbach resonance 
fl4| has also been studied. The work herein is a follow- 
up study of |l3j |. one that includes dissociative damp- 
ing of quantum correlations. In particular, we consider 
a macroscopic atom-molecule superposition created in a 
two step process: first, a joint atom-molecule conden- 
sate is created from an initial atomic condensate with 
strong photoassociation (compared to the s-wave colli- 
sional interaction); second, the photoassociation coupling 
is reduced, and a macroscopic atom-molecule superposi- 
tion is created via strong s-wave collisions; observation 
is based on a Rabi-like oscillation of the system between 
the initial joint atom-molecule condensate and the cor- 
responding macroscopic superposition. Here the main 
concern is to determine whether or not dissociative de- 
coherence prohibits the photoassociative formation of a 
macroscopic atom-molecule superposition. 
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Photoassociation occurs when two atoms absorb a pho- 
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ton, thereby jumping from the two-atom continuum to a 
bound molecular state EE El, 113, E3, E3, [23, 123, 

HE HZl H3- For initially quantum-degenerate atoms, 
there is an analogy in nonlinear optics: photoassocia- 
tion is formally identical to second-harmonic generated 
photons [22|], while s-wave collisions between the par- 
ticles correspond to amplitude dispersion j2^; in each 
case, the inherent non-linearity entangles the particles, 
and potentially leads to macroscopic superpositions. In 
photoassociation, however, the molecular state is excited 
and may spontaneously decay to molecular levels out- 
side the system; also, the molecule may dissociate cither 
back to the atomic condensate, or to a pair of nonconden- 
sate atoms of equal and opposite momentum via rogue 
GlEIHEjl, i.e., unwanted (HiHH, photodisso- 
ciation. These problems can be overcome by coupling 
the excited molecular state with a stable molecular state 
using a second laser. In this two-color photoassociation, 
the population loss via the excited molecular state can be 
made small by applying suitably large intermediate de- 
tunings. Nevertheless, decay of the quantum coherence 
is still possible, even if the populations are unaffected, 
and it is this possibility that warrents an invstigation of 
rogue dissociation in macroscopic atom-molecule super- 
positions. 

The possible creation of a macroscopic atom-molecule 
superposition using photoassociation was reported in 
Calsamiglia et al. Although lacking decoherence ef- 
fects, the results were promising and, moreover, Ref. |l3| 
is an interesting starting point for our study for several 
other reasons: the possibility of a superposition of atom- 
molecule condensates is tempting in itself, since atoms 
and molecules are different objects; the whole dynamical 
process of creating the superposition state can be mod- 
elled; the dissociation environment can be modelled in 
a straightforward manner. In Ref. |14| the dynamical 
approach has been used but, instead of photoassociation, 
the macroscopic superposition state was created with the 
help of a magnetic field and Feshbach resonance. The 
main source of decoherence was the laser-field induced 
interaction between atoms and the electromagnetic vac- 
uum; however, the role of decoherence due to the cou- 
pling between magnetic field modes and the condensate 
is unclear. Also, the study of Ref. E3 considers the case 
of superposition very loosely. One cannot prove that the 
state is a superposition state by using probability dis- 
tributions only. Of course, the possibility of a Rabi-likc 
oscillations to and from the macroscopic superposition 
|l3j. on which observation is based, is more likely in the 
presence of quantum coherence; but, close revivals may 
also occur for nearly statistical mixtures. In all rigor, a 
detailed study of superposition states should include the 
off-diagonal elements of the density matrix, including the 
possibility of dissociative decoherence. 

Hence, we focus on decoherence due to the second or- 
der term of the master equation for the condensate modes 
in the perturbation theory (anomalous quantum correla- 
tions induced by the first order term have been discussed 



elsewhere 29]). Our results are roughly summarized as 
follows. In general, the decoherence timescale tj, should 
be longer than the particular interaction timescale t% , and 
thus macroscopic superpositions are possible only if the 
condition 

& = - » l (i) 

n 

is fulfilled. Our scheme is divided into two phases in 
which different interactions are dominant. In phase J, 
we use photoassociation that is strong compared to colli- 
sions to drive the initial atomic condensate into a partic- 
ular joint atom-molecule which serves as an initial state 
for the superposition engineering. Thus, the time scale n 
is the photoassociation timescale. In phase II, the inten- 
sity of photoassociation laser is turned down, and dom- 
inant collisions between condensate particles drive the 
system into a macroscopic atom-molecule superposition. 
The interaction timescale is then the collision timescale. 
With the given simulation parameters we can express the 
coherence condition as 



sa p[m 3 ](2n 2 + 2n + l)' v ' 

where Pj = 2.712057 x 10 10 , P H = 2.408246 x 10 10 , p 
is the density and n the number of particles in the ther- 
mal cloud. The interesting fact is that the condition 
does not depend on the number of particles in the initial 
condensate, N: With the chosen theoretical setup and 
parameter values, the photoassociation interaction, the 
collision interaction, and the coupling between molecular 
condensate and noncondensate modes depend similarly 
on N; thus, we find that simulations for N = 1000 con- 
densate are applicable to more realistic N = 10 8 conden- 
sates. With values p~2x 10 19 m~ 3 and T ~ 10" 9 K 
we have ^ ~ 3, and the macroscopic atom- molecule su- 
perposition will experience a considerable decay of off- 
diagonal elements. The best chance of avoiding rogue 
decoherence when creating a macroscopic atom-molecule 
superposition is for p — 2 x 10 18 m~ 3 and T ~ 0.1 nK, 
which gives £d ~ 10, but which is just out of reach of 
present ultracold technology |30j . 

The paper is outlined as follows. In Section [H] we 
sketch the photoassociation model given in Ref. [RJ and 
the interactions with the uncorrelated thermal cloud (the 
environment). We calculate the time evolution of the re- 
duced density matrix, i.e., the master equation for the 
condensate modes only. Then, in Sec. IIIII we assign ex- 
plicit values to the model parameters and consider the 
structure of our numerical analysis. The main results 
are also presented there. Section llVI is for summary and 
discussion. 



II. GENERAL THEORY 

Consider a photoassociation laser that removes zero- 
momentum (k = 0) atoms from an initial atomic con- 
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densate |0) and creates molecules in the excited molec- 
ular condensate |1). A second laser couples this excited 
molecular condensate with a stable molecular condensate 
1 2). The annihilation operators of the atomic condensate, 
the excited molecular state and the stable molecular state 
are denoted by ao = a, b and g. Transitions to non- 
condensate atoms arise because a molecule in |1) need 
not dissociate back to the k = atomic condensate |0), 
but may just as well wind up as two atoms of equal and 
opposite momentum (±k) in the state |K) = |k, — k), 
since only relative momentum is conserved. These non- 
condensate atoms are denoted by a k , with ao = a be- 
ing the condensate. The interactions that cause atom- 
molecule and molecule-molecule transitions are charac- 
terized by their respective Rabi frequencies ilk = Oi/k 
and f^2, where /k expresses the wavevector dependence 
of the atom-molecule coupling. The one- and two-photon 
detunings are Ao and 5 = So — ijT, where the sponta- 
neous decay rate, T, of the excited molecular state is in- 
cluded in the one-photon (intermediate) detuning. This 
scheme is illustrated in Fig. ^ We also consider s-wave 
collisions between the atomic and stable-molecular con- 
densates, with interaction strengths A aa , \ gg , and \ ag - 
Anticipating large intermediate detunings, the excited- 
state molecular fraction is low, and its collisions can be 
ignored; similarly, we also neglect all collisions with non- 
condensate atoms. Thus, in the rotating wave approxi- 
mation, the Hamiltonain for the above system is written 
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FIG. 1: Few- level illustration of free-bound-bound photoas- 
sociation. Initially, TV atoms are in the Bose-condensed state 
jO). The free-bound laser then removes two atoms from this 
state and creates an electronically-excited molecule in the 
state |1). The bound-bound laser, in turn, removes an excited 
molecule from the state |1), and creates an electronically- 
stable molecule in state |2). The quasicontinuum of N qc non- 
condensate disociation modes is also shown, where a pair of 
atoms with momentum ±?ik and energy htk = h 2 k 2 /2m are 
taken to occupy the state |K) = |k, — k). The free-bound 
and bound-bound Rabi frequencies are Qk = fii/k and Q2, 
where fk denotes the wavevector (energy) dependence of the 
free-bound coupling. Lastly, the loss rate of the electronically- 
excited molecular state is F, and the respective intermediate 
and two-photon detunings are So and Ao. 



as 

— = -A g^g + Stfb + ^ e k ala k 

k 

_ o XI (^i/kkWa-k + h.c.) - (£l 2 g^b + h.c.) 

k 

+2A ag a t ag t .g + Xggg^g^gg + \ aa ala^aa. (3) 

For simplicity, we take f = fo = F % = ■, 
and /k = /£. We also assume free- particles, so that 
e k = hk 2 /m, where m is the atomic mass. The atomic 
condensate modes ao = a are written explicitly in a mo- 
ment. 

We now develop an effective description by adiabat- 
ically eliminating the excited molecular state from the 
Hamiltonian of Eq. based on the Heisenberg equa- 
tions of motion 

ia = -\fl 1 \e i<t ' 1 a' ! <b + 2\ ag ag^g + 2\ aa a^aa, (4) 
ia k = e k a k --\n 1 \e i ^f k al k b, (5) 

ib = «-^|Oi|e*^/ u aka-k-|0 2 |e-* ff , (6) 

k 

ig = -A g- \fl 2 \e l,p2 b + 2X ag a ji ag + 2X gg g^ gg. 

(7) 

The trick now is assuming that S is the largest frequency 
in the problem, and we can adiabatically eliminate the 
excited molecular state by using b/S ~ 0. Thus, 

= — 25 — 2^ /k"ka-k H ~ 6 g- (8) 

k 

Inserting this into our three-level Hamiltonian of Eq. © 
produces an effectively two-level Hamiltonian. The equa- 
tions can be simplified by denoting the relative phase 
of the lasers with cf> — <f> 2 — fa, and by writing \ — 
|Qi I |fi 2 \/S. This yields 

-r- = -AVff + Xl(^ A o + £ ii)4 a k 

k 

"f E V/kflkfl-k + h.c.) 

k 

+2\ ag a)ag^g + \' aa a)a)aa + \ gg g^g^gg, (9) 

where A' = A + |Q 2 | 2 /<5 and X' aa = A aa — |f2i | 2 / 4<5. Note 
that this is exactly the same form of the Hamiltonian as 
for one-color transitions, but the two-photon Rabi fre- 
quency x nas replaced the one-photon Rabi frequency 

n 1 . 

Next, we proceed to derive the master equation |3lll3^| 
for the stable molecular condensate, applying the same 
ideas as in Ref. j2^. Using the interaction picture, the 
problem is solved to second order in perturbation theory. 
Taking the atom-molecule condensate as the system, and 
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the noncondensate modes as the environment, the Hamil- 
tonian of Eq. (O can be written as 



The dynamics is not altered if, for later convenience, we 
add the constant of motion X' aa (N — N 2 ) to H 2 (see also 
Ref. [13|). Thus, we get H s = H + H T , with 



H ° a t 

T = 3 

-2A<?+aV + (A 99 - 4X' aa )g^gg^g, (12) 

= X^ £fca k ak ' ( 13 ) 



where A = A' - 2X' aa + X gg , and A = 2A^ a - X ag . 
We then apply the Born approximation and calculate 
the master equation of our system in the interaction pic- 
ture. Initially the system |n m )|n a ) and the environment 
Tlk^o l'a)k are not correlated, i.e., p to t = Ps ® pe- The 
equation of the motion of the total density matrix in the 
interaction picture is 



After integration, performing the trace over environ- 
ment and neglecting terms higher than second order, we 
get 



Ps = -~Tr e [H int (t),p s (0)® p e 
n 



dtiTr e [H int (t), [HintiU), p s (0) ® p e }] 



Ui(t) + U 2 (t). 



(15) 



(11) In the interaction picture, our Hamiltonian (|1L)|) becomes 

H int = -h (A^T + rU) , (16) 
where we have used the shorthand notation 

-i(e fc -A)i^ ^y-j 



r = /kOko-k' 

k^O 

2 y 



(18) 



Ptot 



[Hint ; Ptot] ■ 



(14) 



Tracing of Eq. I|15|l results in 

U x (t) = i «T(i)> [A\ Ps (0)] + (Tt(f)> [A, Ps (0)}) (19) 
as the first order term, and the second order term is 



ii 2 (t) = - dh «r(i)r(to> [a\ [a\ p m]] + PHmh)) [a, [a\ Ps (o)]] 

Jo 

-(rt(t)rt( tl ))[A [a, Ps (q)}})- 
i 



(r(t)rt( tl )) [a, Pt 



(20) 



The following six coefficients should be calculated: 



h 
h 

h 
h 
k 



The coefficients 



<r(f)>, 

(rt(t))=jr, 



dti<rt(t)r(ti)), 
dti<r(t)rt(ti)). 



/i = / 2 = / 3 



o. 



(21) 
(22) 

(23) 
(24) 
(25) 
(26) 

(27) 



because for the normal thermalized heat bath the corre- 
lation function (a^a^) is zero, though these have been 



studied elsewhere [29J. In this study, the environment is 
a thermal cloud around the condensates, and therefore, 
it is reasonable to assume that it acts like an uncorre- 
cted thermalized heat bath. The nonzero correlations 
have the form of 



( a iL a -k a k' a -k'> = 4wnl5 k ,k> ■ 
Let us calculate, say, I5: 



(28) 



h = / d*i(rt(i)r(ti)> 

Jo 

= / /k/k'47r^4,fc'e 4 
Ja k k' 

= 4tt / dh J2 fl n l e 

Jo u 



(e fc -A)i e -i(e fc /-A)ii 



i(e Jfe -A)(t-ti) 



(29) 



The sum can be converted into an integral by assum- 
ing that the maximum momentum is very large, i.e., in- 
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finite. This approximation will, however, yield conse- 
quences that are not present in the original system (max- 
imum momentum very large but finite), e.g., irreversible 
dynamics. In our case it means that the coherence of the 
system is decreasing monotonically. Here we have as- 
sumed that the particles are in a free space. Of course, if 
one wants to take the trap into account, one should add 
here the density of the states of the trap. As the first 
approximation, we keep the system as simple as possible; 
neglecting the trap can be justified by the fact that the 
time scales of the trap are much longer than the assumed 
decoherence time scale. 

Next we go to the frequency representation: k = 
y/me/h, de = (2h/m)kdk, so that 

3/2t/ rt poo 

h = ^ l dtl l *Vi/ 2 (e)n 2 ( e )e^-*)^>. 

We assume that the Markov approximation is valid, i.e., 
A is large enough and v / e/ 2 (e)n 2 (e) is a slowly varying 
function in the vicinity of e — A, and thus the correlation 
time scale is so short that we can set t — > oo in the 
integral I4. These conditions are met when HA ~ kT. 
This leads to the principal value integral 

lim / e <(«-A)(t-ti) = w S(e- A) + i — — . (31) 
t^aoJ Q A-e 

Thus, 

** = !! S VZ/2(A)n2(A) + ZA *' (32) 



III. SIMULATIONS 

The aim of the present work is to study the impact 
of the rogue decoherence environment. It is important 
to note that it is a matter of contention whether strong 
two-color photoassociation of a Bose-Einstein condensate 
is even possible [3!J. We insist that it is and, more- 
over, that it is possible to target ground-electronic state 
levels that are sufficiently low-lying to allow neglect of 
vibrational-relaxation losses (see also |3^|) Our simula- 
tions therefore consist of two stages: first, a particular 
joint-atom molecule condensate with real probability am- 



where 



(33) 

While technically related to the Lamb shift, the term A* 
is the shift that appears anytime one connects a bound 
state to a continuum. This has been investigated in the- 
ory |3l EU and experiment [H El OS. The shift is 
neglected in this study. Similarly, for I e we get 



J 6 = ^^^/ 2 (A)[n(A) + lf. (34) 



The master equation in the interaction picture is thus 



Pint = -I 5 [A,[Alp]]~I 6 [A\[A,p}] 

= (I 5 + h)(ApAi +AipA-AAip-pA*A). 

(35) 



The corresponding master equation in the Schrodinger 
picture is 



(36) 



plitude is created in the regime where photoassociation 
is much stronger than collisions; second, the photoasso- 
ciation intensity is reduced and the joint atom-molecule 
condensate evolves under now-dominant collisional inter- 
action into a macroscopic atom-molecule superposition. 

Emipircally [13J, we have found that the parameter 
values that are sufficient to create the proper Phase-/ 
joint atom-molecule condensate are: 

A = 0, (37) 
X = lOViVA, (38) 
= tt/2. (39) 



-[H + H^pg] + p int 

A(g^9P " Pg^g) + T^Xe^aVgp - paJalg) + -xe^+aap - pg^aa) 



- ( X gg - 4A aa)(5 t 55 t 5P - PffWs) + 2A(g t a t 5 ap - pg^atga)] 
+ {h + h){A P A^ +A^pA-AA^p- pA^A), 



where Hj is treated as a perturbation. 
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By proper, we mean a joint atom-molecule condensate 
with a real amplitude which is the origin of the rel- 
ative phase (f> = 7r/2, and where half the original atoms 
have been converted to molecules. The duration of the 
strong two-photon-resonant photoassociation pulse is r, 
which can be obtained from the solution of the semi- 
classical approximation for the molecules ^3], N/A — 
(N/2)tanh 2 (VN X T), which is a result borrowed from 
the theory of second- harmonic-generated photons [27| . 
The condition A = means that the system is on a 
Stark-shifted resonance and the Markov approximation 
becomes dubious. On resonance, the interaction strength 
compared to the evolution time scale increases and one 
expects that decoherence would take a more dominant 
role. Thus, we use a large enough two-photon detuning, 
A > kT, to avoid the resonance and justify the Markov 
approximation: A = 0.1NX, for p — 1. 
The second phase is collision-dominant: 



A = V77x = 0.17VA, 

x = o.iVnx, 

= 0. 



(40) 
(41) 
(42) 



where the atom-atom scattering length is a aa — 5.4 nm 
[ifl Elj . the atom-molecule scattering length is a ag = 
—9.346 nm [Zfl l4ll | . the unknown molecule-molecule scat- 
tering length is approximated as a gg = a aa , the mass of 
a 87 Rb atom is to = 1.443 x 10~ 25 kg, and V is the 
quantization volume that is expedient in a partiuclar 
context (e.g., cubic box, spherical cavity, or harmonic 
trap). Als o, by op tically tuning the atom-atom scattering 
length |3l.l3in4al4ll44ll45|. there is the possibility to re- 
duce the number of parameters by setting X gg — 4A^ a = 0, 
which amounts to a' aa = a gg /8. The effects of the loss 
term r = 12 x 2ir MHz should be negligible if \5\ » 
V; therefore, we choose \S\ ~ 10 3 x 2n MHz. For this 
detuning, and barring an anomalously large molecule- 
molecule scattering length, the light-induced scattering- 
length shift should be possible; however, this really only 
a matter of convenience, and is not neccesary. Finally, 
with the given value of A QQ , the remaining collisional in- 
teraction strength is A = (p/N) x 7.67715 x 10~ 17 Hz, 
where [p] = m~ 3 . 



The collision interactions are 

4:Trha aa 



Aaa — 



A„ 



mV ' 
3irha ag 

mV 
2nha gg 

mV 



(43) 
(44) 
(45) 



We use the molecular number basis: p = \m)(m'\, 
where to and m are the molecule numbers, and N is 
the total particle number. Hence, the coupled equations 
to solve are 



= p m ,m'{«A(m - to') + i2\[(N - 2to)to — (N — 2m')m'] 

2 

- Pm,m' (h + h) \ ( m + 1) 

+ Pm-i.m'^xe-' VCJV - 2to + 1)(N - 2m + 2)m 

- p m , m '+ii\xe-^ VW_ - 2toQ(7V - 2m' - l)(r^T T) 
+ p m+ x, m >i\xe i4 VW_ - 2m)(iV - 2m - l)(m +T) 

- p m . m '-ii\xe i *y/{.N - 2m' + 1)(N- 2m' + 2)m> 

2 

+ Pm-\,m'-\{h + h)\\rmrn' 



+ p m +i, m '+i{h+h)\^{m + l){m' + 1) 



{h + h)\m 
m < N/2 
m > 
to' < N/2 
to < N/2 
to' > 

TO, TO' > 

to, to' < N/2 



'} 



(46) 



Let us stress that these equations of motion (|4*6")) behave 
well at the thermodynamical limit V — > 00 because, as 
typical for photoassociation, X ~ 1/VV h+Ie ~ 

V and p m , m > ~ 1/N. 



The common decay profile of coherence in infinite mod- 
els has the form of e _t / Td , which defines the decoherence 
time Td lJlI- Phase I follows this usual behavior, and 



has the off-diagonal decoherence time 



X 2 (h + h) TO 3 /VVAiV(2n 2 + 2n+l)' 



(47) 

where the total particle number is iV and the fraction of 
noncondensate atoms is n — 1 / [exp(hA / kT) — 1]. Since 
photoassociation is dominant by design, this timescale 
must be compared to the photoassociation time scale 
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x) 1 • The decoherence condition is thus 



W = 



2.712057 x 10 



10 



y/p\^](2n 2 + 2n + 1) 



(48) 



where we have used the detuning A = 0.1 TV A. In most 
simulations we have used a small but realistic temper- 
ature T = 10~ 9 K. Thus, the validity of the Markov 
approximation requires A > 131Hz. The nonzero de- 
tuning, compared with zero detuning, in the phase / 
does not contribute much to the moment r when phase 
I ends. The generalized Rabi frequency 0] is Xgen — 
yjN X 2 + A 2 = 7VAvT00 + 0.01. For simplicity, we do 
not apply the time boundary between phases I and II 
given by the semiclassical approximation [i.e., r defined 
by TV/4 — (N/2) tanh 2 (v A /V"x r )]7 but we instead optimize 
the distance between the two superposition peaks, or, the 
"size" of the superposition, to be a large as possible. An 
optimal superposition size, and equally-sized probability 
peaks, are obtained only if the phase I ends at the mo- 
ment r when Tipg^g ~ 0.296TV/2. This end point is very 
sensitive: the accuracy in molecular fraction should be 
~ ±0.001 in order to get equally sized peaks. Also, the 
end point depends a bit on the particle number. The 
above mentioned result is valid for N = 1000. 

Phase II follows similar behavior, having the same 
decoherence time as in Eq. (|47|) . However, since colli- 
sions are dominant in phase //, we compare the deco- 
herence time in phase 77 with the collision time scale 
t c ~ (2TVA) -1 . Hence, the decoherence condition is 



r d 

Sd.IJ — — 

Tr. 



2.408246 x 10 1 



v /p[m 3 ](2n 2 + 2n + 1) 



(49) 



If phase / results in suitable initial conditions to pro- 
duce a macroscopic superposition in phase II, decoher- 
ence in phase II is the only remaining thing to prevent 
the superposition state. The correlations created by colli- 
sions are the mechanism that, in ideal conditions, results 
in a macroscopic superposition of atomic and molecular 
condensates. Therefore, if the decoherence time scale is 
smaller than the collisional time scale then decoherence 
dominates over collisions and prohibits the creation of 
the macroscopic superposition. The constraint to the su- 
perposition here is £d,li >> 1- 

There are several damping scenerios to consider; al- 
though, due to limited computational resources, we con- 
sider the case of TV = 1000 particles. In scenario A there 
is no decoherence present, i.e., £,d,i — £,d,n ~ oo. In 
scenerio B, the parameter values p = 1.88 x 10 19 m -3 
and T = 10 -9 K result in moderate decoherence, i.e., 
£d,i = 3.44 and £d,n = 3.05. The scenario C is the 
borderline scenario with £d,i = £d,n = 1- Correspond- 
ing temperatures and densities can be calculated from 
Eqs. gSJl and gjj. In scenario D, we assume a tem- 
perature T = 10~ 10 K and density p = 1.88 x 10 18 m- 3 , 
he., £,d,i = 10.88 and £d,ll = 9.65. The results are illus- 
trated in Figs. 12141 In Fig. [21 we present the probability 
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FIG. 2: The probability distributions of scenarios A (no 
damping), B (moderate damping), C (strong damping) and D 
(weak damping) at (a) \t\ — 8.6 x 10~ 5 when the stage I ends 
and (b) at At2 = 1.0486 x 10 -2 when the superposition state 
emerges at the end of stage II. Decoherence obviously affects 
the probability distributions: stronger decoherence makes the 
probability peaks lower and wider. Recall that the stage / 
and II decoherence conditions for the various scenerios are 
(A) Ui = Un « oo; (B) Ui = 3-44 and fan = 3.05; 
(C) Ui = Un = 1; and, (D)£ dJ = 10.88 and Un = 9.65. 



distributions of all scenarios at two different instants of 
time. The first instant (ti) is when the phase I ends, and 
a particular joint-atom molecule condensate with a real 
probability amplitude is created. It serves most ideal ini- 
tial conditions to begin the collision dominant phase //. 
The second instant fa) is when the superposition peaks 
should emerge. Figure presents the density matrices of 
all scenarios at ti, and Fig. 21 the density matrices of all 
scenarios at t%. 

Overall, while decoherence is strong enough, the quan- 
tum correlations are damped out and sharp peaks in 
probability distributions smoothcn. The decay of off- 
diagonal elements seems to be considerable also in the 
regime of moderate decoherence 1 < £<j < 6. Strong de- 
coherence effects £d < 1 would be even faster. In order 
to get a macroscopic superposition with the given setup, 
one should be safely in the regime of weak decoherence, 
i.e., £d > 10. However, the regime is not reachable using 
parameter values consistent with present-day technology. 
Of course, there are a few possibilities to improve the re- 
sult and get £^ > 10. First, one can reduce the density 
p, but 10 19 m -3 is somewhat rarified as it is, even for 
a condensate, and achieving 2 x 10 18 m~ 3 could present 
challenges. The second possibility is to enter the non- 
Markovian regime, which could considerably improve the 
t;d at least in phase /; however, the non-Markovian case 
will be studied elsewhere. Other tricks are related to 
manipulating the interaction or the environment. The 
symmetrization of the environment |12| would be a hard, 
maybe impossible, task to apply here due to the inherent 
nonlinearity; but, creating specially a correlated environ- 
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Fig 3a Fig 3b 




FIG. 3: The density matrices of scenarios A (a, no damp- 
ing), B (b, moderate damping), C (c, strong damping) and 
D (d, weak damping) at the end of stage I at ti. Deco- 
herence effects are clearly visible in scenarios B and C, as 
the decay of off-diagonal elements has resulted in ellipsoid- 
gaussian wave packet. Again, recall the decoherence condi- 
tions: A, £ diI = (, d>II ~ oo; B, £ di i = 3.44 and £ d ,n = 3.05; 
C Ui = Uu = 1; and D, = 10-88 and £ d , n = 9.65. 
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FIG. 4: The density matrices of scenarios A (a, no damp- 
ing), B (b, moderate damping), C (c, strong damping) and 
D (d, weak damping) at ti-, when a macroscopic superposi- 
tion state should emerge. Only a small fraction of off-diagonal 
elements have survived from decoherence, even in the sce- 
nario D, that has extreme parameter values of T = 10~ 10 K 
and p = 10 18 m" 3 . A gain, recall the decoherence conditions: 
A, U,i = £d,« ~ oo; B, = 3.44 and ^ d>n = 3.05; C 
id,i = £d,n = 1; and D, £ djI = 10.88 and £ dJI = 9.65. 

merit might help since an uncorrelated heat bath is quite 
a harsh environment for coherence phenomena. 

Finally, we note that the observation of the possible 
superposition state could present a problem. The as- 
sumption that the coherence of the state after phase I 
can be verified, and thus one only needs to consider the 
density matrix profiles of the revival of this initial state 
after the superposition state [ill. Il4| , is perhaps too sim- 



ple, at least in the Markovian regime. While the system 
is driven from superposition to the revived state, rogue 
dissociation is still damping coherence out of the system. 
Also, coherent tricks, such as phase-imprinting the re- 
vived state, would need extra time, and this introduces 
more decoherence. The best candidate for observing the 
macroscopic superposition is probably to try to entangle 
the superposition state with another system, and then to 
apply Bell-type correlation experiments. 



IV. DISCUSSION 

The important point is that our approach is a dynam- 
ical one in which we make the macroscopic superposition 
in the presence of the decoherence, which differs from the 
standard way in which one takes the macroscopic super- 
position for granted and then puts it into a decoherent 
environment and looks how long it would last there. Such 
studies do not make claims whether it is possible to cre- 
ate a macroscopic superposition or not, but they only 
make claims whether the initially prepared macroscopic 
superposition can be observed or not. Therefore, our 
photoassociation-based approach suggests a possible ex- 
perimental creation of the macroscopic superposition in 
a decohering environment. 

Hence, the aim of this study is to understand the ef- 
fects of the rogue-dissociation-induced decoherence while 
trying to create a macroscopic superposition of Bose- 
condensed atoms and molecules. In particular, ours is 
a two-stage scheme: in the first stage, the photoassocia- 
tion coupling is strong compared to the s-wave collisional 
coupling, and a specific joint atom-molecule condensate 
is created; in the second stage, the photoassociation cou- 
pling is turned down such that it is weak compared to 
the collisional coupling, and the system evolves into a 
macroscopic atom-molecule superposition. There is a 
constraint ^ >> 1, which should be satisfied in order to 
get a macroscopic superposition of atoms and molecules. 
The dependency on relevant system parameters of the 
constraint ^ was calculated in both stages of our theo- 
retical setup. An interesting fact is that, for the given 
parameter ratios, £<j does not depend on the number of 
particles in condensate N. Results for a modest, say, 
N = 1000 particle superposition state will therefore ap- 
ply to a N — 10 9 particle superposition state with the 
same density, given properly scaled laser frequencies and 
intensities. Using ^ one can easily, without heavy simu- 
lations, consider whether the macroscopic superposition 
state of atoms and molecules is within reach or not with 
given experimental setup parameters. 

Without any tricks, the best chance to defeat rogue 
decoherence in creating an atom-molecule superposition 
corresponds to a density p = 2 x 10 18 m -3 and temper- 
ature T ~ 0.1 nK. Nevethless, while this temperature 
may be accessible by making the trap shallow, the off- 
diagonal elements of density matrix will still experience 
considerable decay. One way to get better macroscopic 
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superpositions in the present scheme is to challenge the 
assumptions or approximations of our study or, in other 
words, engineer the environment. We assumed that the 
environment is an uncorrelated heat bath, the system 
parameter values should be in the regime where Markov 
approximation is valid, the environment particles are free 
and the environment is infinite. Indeed, it may well prove 
interesting to examine the non-Markovian regime, where 
a greater part of the physically-relevant parameter space 
is available, and/or to use a specially correlated environ- 
ment. 

Rogue dissociation is an important source of decohcr- 
ence that should not be neglected in coherence studies 
of atom-molecule Bose-Einstein condensates. Unfortu- 



nately, it seems to push coherent atom-molecule engineer- 
ing (macroscopic superpositions, quantum computation, 
etc.) one step further away from realization. 
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